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CO , We compute the full bispectra, including cross bispectra, of primordial curvature and tensor per- 

turbations in the most general single-field inflation model whose scalar and gravitational equations 
'" of motion are of second order. The formulae in the limits of k-inflation and potential-driven inflation 

\mS \ are also given. These expressions are useful for estimating the full bispectra of temperature and 

■ polarization anisotropies of the cosmic microwave background radiation. 
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$H 1 I. INTRODUCTION 

i ^ , ' The non-Gaussianities of the temperature and polarization anisotropies of the cosmic microwave background (CMB) 
radiation now receive increasing attentions because they are important tools to discriminate models of inflation [l|, Q ■ 
t— I ■ Ongoing and near future project such as Planck satellite , CMBpol mission Q , LiteBIRD satellite @ would reveal 
the properties of the temperature and polarization anisotropies in detail. Such E-mode polarization anisotropies are 
sourced by both curvature and tensor perturbations Q, while only tensor (and vector) perturbations can generate 
B-mode polarization anisotropies 0. 1 Therefore, even when one estimates the "auto" bispectra of the temperature 

■ and the E-mode polarization fluctuations, not only the auto bispectra but also the cross bispectra of the primordial 
curvature and tensor perturbations are indispensable. 

For slow-roll inflation models with the canonical kinetic term Q, Maldacena evaluated the full bispectra, includ- 
' ing the cross bispectra, of the primordial curvature and tensor perturbations [p. Inflation models are now widely 
generalized into more varieties such as k-inflation DBI inflation 11 1, ghost inflation[l2j, G-inflation [l3|], and so 
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on. However, almost all the works on the non-Gaussianities in these inflation models concentrate only on the auto 
, | bispectrum of the curvature perturbations p^4l6| , which is insufficient for evaluating the bispectra of the temperature 
■ and E-mode polarization anisotropies of the CMB, as explained above. To our surprise, as far as we know, the full 
I bispectra of the primordial curvature and tensor perturbations have not yet been obtained even for k-inflation [Toj 
except for Ref. [l?} where the primordial scalar-scalar-tensor cross bispectrum has been calculated for inflation models 
with an arbitrary kinetic term. There are several related works on the primordial cross bispectra. In Ref. [lH , the au- 
thors show the primordial tensor-scalar cross bispectra induced from a holographic model and the scalar-scalar-tensor 
correlation has been discussed in the calculation of the trispectrum of the scalar fluctuations fl9j , so-called "graviton 
exchange" , and also in the context of one-loop effects of the scalar power spectrum [2(| . In Ref. [2l[ , the authors 
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1 Though vector perturbations can also generate both E-mode and B-mode polarization anisotropies, they only have a decaying mode in 
linear theory and hence suppressed in the standard inflationary cosmology based on scalar fields. 
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calculate the correlation between primordial scalar and vector (magnetic fields) fluctuations in possible inflationary 
models of generating primordial magnetic fields. 

Among the inflation zoo, the generalized G-inflation model [22| occupies the unique position in that it includes 
practically all the known well-behaved single inflation models since it is based on the most general single field scalar- 
tensor Lagrangian with the second order equation of motion, which was proposed by Horndeski more than thirty years 
ago [23| and was recently rediscovered in the context of the generalized Galileon 0, HEf . Indeed, it includes standard 



canonical inflation [I], non- minimally coupled inflation [2q ] including the Higg s inflation 27 [.extended inflation 



28], k-inflation [Io|, DBI inflation [ff|, R 2 inflation [IHH, new Higgs inflation |3C§ , G-inflation [I3| , and so on. Thus 



once we analyze properties of the primordial curvature and tensor perturbations in the generalized G-inflation, one 
can apply the result for any specific single-field inflation models. 

So far, the power spectra of scalar and tensor fluctuations were studied in [22j and the general formulae for them 
have been given there. It has been pointed out that the sound velocity squared of the tensor perturbations as well as 
that of the curvature perturbations can deviate from unity. Then the auto bispectrum of the curvature perturbations 
was estimated in Refs. [HI [32[ (see also (33|) and found to be enhanced by the inverse sound velocity squared and 
so on. More recently, the auto bispectrum of the tensor perturbations was investigated in Ref. [34j ] and found to be 
composed of two parts. The first is the universal one similar to that from Einstein gravity and predicts a squeezed 
shape, while the other comes from the presence of the kinetic coupling to the Einstein tensor and predicts an equilateral 
shape. What remains to be studied are the bispectra of the primordial curvature and tensor perturbations in the 
generic theory. 

In such a situation, in this paper, we compute the cross bispectra of the primordial curvature and tensor perturba- 
tions in the generalized G-inflation model. The formulae in the limits of k-inflation and potential-driven inflation are 
also given as specific examples. 

The organization of this paper is given as follows. In the next section, we briefly review the most general single 
field scalar-tensor Lagrangian with the second order equation of motion. In Sec. Ill, quadratic and cubic actions for 
the primordial curvature and tensor perturbations are given. The full bispectra, including the cross ones, for them 
are discussed in the section IV. The special limits for them in the cases of k-inflation and potential driven inflation 
are taken in Sec. V. Final section is devoted to conclusion and discussions. 



II. GENERALIZED G-INFLATION — THE MOST GENERAL SINGLE-FIELD INFLATION MODEL 

The Lagrangian for the generalized G-inflation is the most general one that is composed of the metric and a 
scalar field (f> together with their arbitrary derivatives but still yields the second-order field equations. The Lagrangian 
was first derived by Horndeski in 1974 in four dimensions [23j . and very recentl y it was rediscovered in a modern form 
as the generalized Galileon (24|, i.e., the most general extension of the Galileon j3ol[36|. in arbitrary dimensions. Their 
equivalence in four dimensions has been shown in Ref. [22| . The four-dimensional generalized Galileon is described 
by the Lagrangian: 

C = K(<j > ,X)-G 3 (4>,Xp4, + G 4 {ct>,X)R + G 4 x [(n</>) 2 - (V M V„0) 2 ] 

+G 5 (0, X)G^W^4> - \g 5X [(D0) 3 - 3D0(V M V^) 2 + 2(V Al V„0) 3 ] , (I) 

where K and Gi are arbitrary functions of (j> and its canonical kinetic term X := — (50) 2 /2. We are using the notation 
Gix for dGi/dX . The generalized Galileon can be used as a framework to study the most general single- field inflation 
model. Generalized G-inflation contains novel models, as well as previously known models of single-field inflation such 
as standard canonical inflation, k-inflation, extended inflation, and new Higgs inflation, and even R 2 or f(R) inflation 
(with an appropriate field redefinition). The above Lagrangian can also reproduce the non-minimal coupling to the 
Gauss-Bonnet term 1221]. 



III. GENERAL QUADRATIC AND CUBIC ACTIONS FOR COSMOLOGICAL PERTURBATIONS 

In this section, we present the quadratic and cubic actions for scalar- and tensor-type cosmological perturbations 
based on the most general single-field inflation model. Employing the Arnowitt-Descr-Misner formalism, we write the 
metric as 



ds 2 = -A 2 dt 2 + g {j (dx l + N l dt) (dx j + N j dt) , 



(2) 
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where 

N=l + a, Ni = dip, gij = a 2 (t)e^(e h ).., (3) 

and {e h )ij = <% + + (l/2)hikhkj H — • • We work in the gauge in which the fluctuation of the scalar field vanishes, 
= <t>(t). Concerning the perturbations of the lapse function and shift vector, a and (3, it is sufficient to consider the 
first order quantities to compute the cubic actions, as pointed out in The first order vector perturbations may be 
dropped. The curvature perturbation in generalized G-inflation is shown to be conserved on large scales at non-linear 
order in [37j |. 

Substituting the above metric to the action and expanding it to third order, we obtain the action for the cosmological 
perturbations, which will be written, with trivial notations, as 

S = J (C hh + C ss + C hhh + C shh + C ssh + C sss ) . (4) 

The first two Lagrangians are quadratic in the metric perturbations, which have already been obtained in Ref. [22j. To 
define some notations used in this paper, we will begin with summarizing the quadratic results in the next subsection. 
The third and last cubic Lagrangians have been derived in Refs. [34| and [HI, [Hf, respectively, but for completeness 
they are also replicated in this section. The mixture of the scalar and tensor perturbations, C s hh and C ss h, are 
computed for the first time in this paper. 



A. Quadratic Lagrangians and primordial power spectra 

The quadratic terms are obtained as follows [j^ . 



1. Tensor perturbations 



The most general quadratic Lagrangian for tensor perturbations is given by 



£hh = — 



• 2 J~T 



where 



•Ft 
Gt 



Gi — X i^pG^x + G^cjy 
G 4 — 2XG^x — X I H(f>G 5 x — G 5ctl 



(5) 

(G) 
(7) 



Here, a dot indicates a derivative with respect to t, Gi$ := dGi/d<j) and the propagation speed of gravitational waves 
is defined as := Ft/Gt- The linear equation of motion derived from the Lagrangian ([S]) is 

E% := d t (a 3 GTh l0 ) - aF T d 2 h l0 = 0. (8) 

In deriving the above equations, we have not assumed that the background evolution is close to de Sitter. They can 
therefore be used for an arbitrary cosmological background. 
We now move to the Fourier space to solve this equation: 



hij(t,x) 



d 3 k 



hij ;(t, k)e ik x . 



(9) 



It is convenient to use the conformal time coordinate defined by drj ~ dt/ 'a. We approximate the inflationary regime 
by the dc Sitter spacetime and take Tt and Gt to be constant. The quantized tensor perturbation is written as 



Kr ? ,k)=^[/ lk (r ? )eg ) (k)a s (k) + / l *_ k ( 7? )e*f ) (-k)at(-k) 



where under these approximations the normalized mode is given by 

iy/2H 



\jT T c h k 3 



(l + ic h k V ) e~ lChk7 K 



(10) 



(11) 



Here, is the polarization tensor with the helicity states s = ±2, satisfying e\f (k) = = kje\j(k). We adopt the 
normalization such that 



eW(k)e*f s '>(k)=^, 
and choose the phase so that the following relations hold. 

e^ ) (k) = e^ ) (k) = eg ) (-k). 
The commutation relation for the creation and annihilation operators is 

[a s (k),al,(k')] = (2 7 r) 3 ( 5^ ( 5(k-k'). 
The two-point function can be written as 



where 



{h tJ (k)h kl (k')) = (2n) 3 S(k + k')T tjM (k), 
P« lW (k) = \h k \ 2 U ijtk i(k), 



n i , )W (k) = ^eg ) (k) e ^ ) (k). 



(12) 



(13) 



(14) 



(15) 
(16) 



(17) 



The power spectrum, Vh = (k 3 /2ir 2 )'Pij,ij , is thus computed as 



2 H 2 



IT 2 T T c h 



Chkri= — 1 



(18) 



2. Scalar perturbations 



The quadratic Lagrangian for the scalar perturbations is given by 



where 



e 



C ss = a 3 



XK X + 2X 2 K XX + 12H<j>XG 3X + 6H<j>X 2 G 3X x - 2XG 34> - 2X 2 G 



6H 2 G 4 



+6 



H 2 (JXG4X + l6X 2 G 4 xx + ^X 3 G 4 xxx) — Hip + 5XG 4( px + 2X 2 G 4( f > xx) 

5c/>X 



-30H 3 ^XG 5X + 26H 3 <PX 2 G 5XX + ±H 3 ^X 3 G 5X xx - 6H 2 X (6G 50 + 9XG 5 ^ X + 2X 2 G b<t>X x] 
-4>XG 3X + 2HG 4 - 8HXG 4X - 8HX 2 G 4X x + 4>G 4c/l + 2X<j)G 44>x 



-H 2 j> (5XG 5X + 2X 2 G 5X x) + 2HX (3G 50 + 2XG 5<j>x ) • 
Varying Eq. (|TT))) with respect to a and /3, we get the first-order constraint equations: 

Ea - ^d 2 p + 36C - %d 2 C = 0, 
Oa-g T ( = 0, 



(19) 



(20) 
(21) 

(22) 
(23) 



which arc solved to yield 
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with V := d~ 2 (. Plugging Eqs. ([24]) and ([HJ to Eq. ([19]), wc obtain 



C ss = a? 



' 2 **S 

SsC — 5- 00 



where wc have defined 



•Fs '■= -—[~Qt)^^t, 



1 d / a 
adi ve 1 

E 2 

ys := Q2 + 3£t- 



(26) 

(27) 
(28) 



The sound speed is given by c 2 := J^s/Gs- The linear equation of motion derived from the Lagrangian (|26ll is 

E s := d t (a 3 GsC) - aT s d 2 C = 0. (29) 

The scalar two-point function can be calculated in a way similar to the case of the tensor perturbations. We move 
to the Fourier space: 

d 3 fc 



C(*,x) 



(2tt): 



■C(i,k)e 



ikx 



(30) 



and proceed in the de Sitter approximation, assuming that Fg and Qs are almost constant. The quantized curvature 
perturbation is written as 

cfa,k) = Ck(»?)o(k) + e k w« f (-k), 

where the normalized mode is given by 

m 



{1 + icsh]) e- c ° k11 . 



2^F s c s k A 

The commutation relation for the creation and annihilation operators is 

[a(k)y(k')] - (2^) 3 <5(k-k'). 
Thus, the power spectrum is calculated as 

(C(k)C(k')) = (2^) 3 5(k + k')|fn 



1 H 2 

V< ~ %n 2 F s c s 

From Eqs. (fTS)) and (|3"5"]) . tensor-to-scalar ratio r is given by 



c s krj— — l 



16 



(31) 

(32) 

(33) 

(34) 
(35) 

(36) 



where we have assumed that the relevant quantities remain practically constant between the horizon crossings of 
tensor and scalar perturbations that occur at different time in case Ch ^ c s [38j . 

B. Cubic Lagrangians 

We now present the most general cubic Lagrangians composed of the tensor and scalar perturbations. Wc would 
like to emphasize that in deriving the following Lagrangians the slow-roll approximation is not used. 



1. Three tensors 



The Lagrangian involving three tensors was derived in Ref. [34{: 



Chhh = a 



where we defined 



/i • • J~t ( 1 
—//,,//,/,///,, + — - I hikhji — -hijhki ) hij^ki 



/j, := <j)XG 5X - 



(37) 
(38) 
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2. Two tensors and one scalar 



The interactions involving two tensors and one scalar are given by 



■'shh 



3g T 
8 

Gt 



^Ch 2 - 



Qt 



Ji 

2a 



1 ■ . 



h:h: 



where 



T := 2G 4 - 8XG iX - 8X 2 G iX x - 2H4> (5XG 5X + 2X 2 G 5 xx) + 2X (3G 54> + 2XG 5< px) . 



This quantity can also be expressed in a compact form T = dO/dH. 

Substituting the first-order constraint equations to Eq. (|39|). the Lagrangian reduces to 



C-shh — a 
where 



biChfj + -^(hij t khij t k + l >::i.-.i.'h,J>,,j, + b^h 2 ^ + — |c? 2 C/ifj + b & ^^h ik h jk + -^(,ijh ik h jk 



In 

b 3 
In 
h 
h 
h 



ml Qt d ( Gt \ 
9J T 3 df \OTtJ 



3gr 
8 

Fs_ 

8 ' 

4 ' 
Gt 
86Jt 

fJ-Gr f J^sGt 

le" \TtGs 

2Gt 
V Gt 

2 e ' 



(G 2 T TF T ) + ? 



Gs_ 

Gt 



=£- - 1 



HQ 2 



Gs 
HQs 



G$d_ 

4 dt \QF T 



and 



E. 



hh 



H G 



T i 



hiE s 



4Gs O^t lJ ' 20T T \2 ' g T 
The last term E s hh can be removed by redefining the fields as 

2 



C 



c c + 



8fe e j- T 



hh - 



(39) 
(40) 

(41) 

(42) 

(43) 
(44) 

(45) 

(46) 
(47) 
(48) 

(49) 

(50) 
(51) 



The contribution to the correlation function is however negligible because the above field redefinitions involve at least 
one time derivative of the metric perturbation, which vanishes on super-horizon scales. 



3. Two scalars and one tensor 



The interactions involving one tensor and two scalars are given by 



3Gt 



20a '' ,,//,, + -a + —t" i, j I >,,,.!,!, - 2g T (P,ijhij + fj£P,ijhij - FrCiQh 



(52) 
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Substituting the constraint equations, we obtain the reduced Lagrangian: 



£-ssh = a ' 

where 

ci = 

C2 = 



Cl 



C2 



°5 o2, 



ss/i i 



and 



with 



— (T S -T T ) + — 



1 i?r 

2 + 46 



Qt 
HQt~ 



Id 
4dt V 9 



M-^S , 2HQ T ^ d 



Qt 



(-) 



C4 



Co 



Cfi 



r 



— i - 



2 dt V 29 St/ 



3H 



Qt 
Qt 



V_ + JL 
26 Q T 



Q\ - TT T 2Hfi 



2eg T 



Q T \Q\- TT t 2Hh d 



26 
Ql_ 



20Q T 
<7t 



6 



6 dt ve 
(I) 



dt 



^ (3J- T - J*) 



dt 



IQt^t, ( ^ 



E ss h — /<9 2 d%E a + fijEjj, 



r 

26 

9<?T V 2 £?T 



g 2 tij&kij, 



Qt 



The field redefinition: 



(53) 

(54) 
(55) 

(56) 
(57) 
(58) 
(59) 

(60) 

(61) 
(62) 

(63) 
(64) 



removes the last term E ss h- Since all the terms involve at least one derivative of the metric perturbation, the field 
redefinition docs not contribute to the correlation function on super-horizon scales. 



4- Three scalar s 

For completeness, here we give the cubic Lagrangian for the scalar perturbations derived in Refs. [H], [Hj]. The 
cubic Lagrangian for the scalar perturbations is given by 



r — 



where 



3SC + SC + (F - Qt) % - 



-— (E + 2XT,x + HE) a 3 + a 3 
o 

+AafiaCC,ii ~ (A«A« _ AwAj'i) + ~ a (AtfCw 



a 2 - 2aQaQP >i + 18a 3 6aCC 



aQ T aOO + 3a d TaC z + 2a A nC + aTKC.iC.i " ^Q T CC + 2aQ T P^ - 2a^ u C + 2aQ T p M & 
1 /3 



+ -1^5 -2%A* 



(65) 



l2(/)XG 3 x + 6(t>X 2 G 3X x - \2HG A 



-6 



2H (7XG 4X + l6X 2 G 4X x + 4X 3 G 4X xx) - <f> (G A<j> + hXG i<t>x + 2X 2 G 



i<j,XX J 



90H 2 (j)XG 5 x + 78H 2 <j)X 2 G 5 xx + \2H 2 ^X 3 G bX xx - 12HX (6G 5l p + 9XG 5 ^ X + 2X 2 G^ X x) ■ (66) 
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Using the first-order constraint equations to remove a and (3 from the above Lagrangian, we obtain the following 
reduced expression: 



C sss = / dtd A xa A Gs 



6H 



C + C2CC + C 3 ^-C (diO + 2C 4 CftC3V + 2C 5 d 2 C (dupy 



(67) 



with ip = d 2 £. There are five independent cubic terms with coefficients: 



Ci = - 



2H 2 



2~G T 3Q% 



39 3 £ s QT S 
2H 

, 3Gt I Gs , 3Gt 



r {Gs - 2T S ) + 36^ {Gt - Gs) + ^Gt (2Gt - Gs) 



9 2 STs 



9T 

e 



f f \ 2(S-XS x )g£ EGt {3Gt 1 
&tl[ GS-G- T ) + WGS + e^U"" 1 



e V F s Gs 



- 1 



3^ 



6H 3 GsGt ( r ?Gt\ 



(68) 



6 2GT S J 



3/i + ^e h 



(69) 



H 



(3Gs - %Gt) Gt fJ-Gs FGs 



49Js 



2Gl 



AQGt 



H 



Gs ( G T Tg T \ H 3 GsG. 



TGt 
29 



(70) 



Gt \ H 2 Gs 



TGt\ 



Wt + 3HGs {& + iek - 2#> ) + 3 e^f ^ + W ) 



(71) 



3£s 3iJ£s / fx 



8G T <±Gt \Gt 26 



(72) 



IV. PRIMORDIAL BISPECTRA 

Having obtained the general cubic Lagrangians composed of the scalar and tensor perturbations, we now compute 
the bispectra in this section. Here, we use the mode functions in exact de Sitter. 



A. Three tensors 



Let us consider three-point function of the tensor perturbations: 



} (hhh) 



(h idl (ki)k aia (k 2 )h ish (kg)) - (27T) d J(k! + k 2 + ks^V 

_ (2^M^(ncw) 



i1j1i2.72i3.73 



A 



(GR) 



1,31,31,3 V l lJl l l3ll3j3 ' iljli2j2i3j3 y ' 
ft- l"'2"'3 



(73) 
(74) 



where -4, ™ .■ ,■ „■ and A- ■ / ■ - ■ represent the contributions from the h 3 term and the h 2 d 2 h terms, respectively. 

'Ul'2j2«3j3 *1J1*2J2«3J3 ^ 



f (now) and ^. (GR) 
Each contribution is given by 

H jjL k-y k^ k^ 



■^iiniinizh ~ J(3 ^i 1 j 1 ,lm(kl)Ili 2 j 2im n{k2)IIi 3 j s ,nlQ i 3), 



r(new) 

•^■ifhlhish A { "< ' >ik ( k l ) n »2 j 2 ,ii ( k 2 ) 
where i-C = fox + fc 2 + ^3 and 



fc3fcfc3ini3j 3i ij(k 3 ) — -k 3 ik3kHi 3 j 3 ,jl(k 3 ) 



A(k 1 ,k 2 ,k 3 ) := - — 
16 



fcifc 2 fc 3 



(75) 
, (76) 

(77) 
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The first term •^ < i^™i 2 j 2 i 3 j 3 is proportional to G$x and hence vanishes in the case of Einstein gravity, while the second 

term ■Ai 1 j 1 i 2 j 2 i 3 j 3 is universal in the sense that it is independent of any model parameters and remains the same even 
in non-Einstein gravity. 

In order to quantify the magnitude of the bispectrum, we define two polarization modes as 

fW(k):=fty(k) e ;W(k), (78) 
and their relevant amplitudes of the bispectra as 

(^(k^Hk^Hk,)) = (2nf5(k 1 +k 2 + k 3 )^3 +^(GR) 3 ) • (79) 

From Eqs. ([75)1 and ([To]) , the amplitudes *4(new) 3 (GR) are easn y calculated as 34 1 

= ^^F(s^,s 2 k 2 ,s 3 k 3 ), (80) 

-%r) 3 = y (sih + s 2 k 2 + s 3 k 3 f Fisth, s 2 k 2 , s 3 k 3 ), (81) 



where 



F(x,y,z) := 77 2 \ 2 ( x + V + ^) 3 (z ~V + z ){ x + V~ z )(x -y-z). (82) 
o4 ar?/ z z z 



As pointed out in Ref. [34J, has a peak in the equilateral limit, while A^q^) m the squeezed limit. 
It would be convenient to introduce nonlinearity parameters defined as 

TS1S2S3 

fSlS2S 3 _ Qn (now),(GR)fcl=fc2=fc3 foo\ 

^NL(new),(GR) ~ JU K 3 ' ^ 00 <' 

which arc quantities analogous to the standard /nl for the curvature perturbation. We find 

/NLTnow) = "15^68 [3 + 2(SlS2 + 5253 + S3Sl)] ^' (M) 



or, more concretely, 



f +++ = T++- = 5 _ElL (SK\ 

•'NL(new) n52 Q T ' •'NL(new) 10368 g T ' > 

with /NL(ncw) = /nl(»cw) and /NL(rTcw) = /nl(L) ■ ( This symmetry arises because parity is not violated.) As for 

■S1S2S3 

NL(GR 



/nugrp we have 



so that 



~ 85 

/nl ( gr) = [21 + 20( SlS2 + s 2 s 3 + sa*!)] , (86) 

J+++ - J - 255 J++- - 7+ 85 (o 7 \ 

•'NL(GR) ~~ J NL(GR,) ~~ 1024' ■'NL(GR) ~~ ■> NL(GR) ~~ 27648' 

As defined in Eq. ([74]), B^] jai3j3 is normalized by V\. This normalization can be justified when one concentrates 
on the non-Gaussianity of the B-mode polarization. Because the B-mode polarization can be generated by not 
curvature perturbations but tensor perturbations (except for lensing contribution), the size of the non-Gaussianity of 
the B-mode polarization could be directly characterized by /NL S ( 2 n s c 3 w ) (qr) • 

However, it should be noticed that tensor perturbations can generate not only the B-mode polarization but also 
the temperature fluctuation and the E-modc polarization. The latter two are mainly generated by the curvature 
perturbations. Therefore, when one would like to quantify the auto and cross bispectra of the temperature fluctuation 
and the E-mode polarization, it would be better to normalize Biifcilj 3 i s j s by V^, namely, 



T3pp" 

.(new),(GR) _ 2 T(ncw),(GR) , . + „ + T,. c„w^ „ m „ /isis 2 s 3 



R (hhh) _ I C 

iijiiahhh ~ )L3L3jL3 
"' 2 "'3 



/ . (new) a(GK) \ /no\ 



where A\ in ^ 2l J 3 = r 2 ^-^^ jvith r being the tensor-to-scalar ratio. In the same way, -4 (now)i(GR) 

Jl /)SiS 2 «3 anr \ fSlS 2 S 3 _„2fSlS2S3 

(new), (GR) u •> NL( new) , (GR) •> NL( new),(GR) ' 
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B. Two tensors and one scalar 



The cross bispectrum of two tensors and one scalar is given by 

J_1,J!, 



(C(ki)Mk 2 )Mk3)> = (2^ 3 ^(ki +k 2 + k 3 )B i ( 7 C ^ ) (k 1 ,k 2 ,k 3 ), (89) 



where -By ^ is of the form: 



Bffi = , 3 1 3 x- S- , E ^^(ki.ka, ka)l(') (fcx, fca, fc 3 ) + (k 2j i, j o k 3 , fc, /). (90) 



Each contribution is given by 



;(D — tt. . n,_YTT.. n,^ _ u_ . vO) _ ki -k 3 ^(i) 



Vij,ki ~ n ij, mn (k 2 )n feii „ m (k 3 ), V ijM - k 2 • k 3 V ij fei , V y fci - — ^ — V ij,ki 



and 



vgL=vW ;) Vg^^vgi,, Vy 6 i ; = fc lm fc ln n ijW (k 2 )n fe , nm ,(k 3 ), vgy = A^g,, (91) 



(1) = 1 cjkjkKch + A") 
ff 2 A-' 2 



I 1 



(2 ) 1 c^fcf + 2c 2 s c h k 2 (k 2 + fc 3 ) + 2c s c 2 h k 1 (k 2 + fc 2 fc 3 + fcf) + 4(fc 2 + fc 3 )(fc 2 + fc 2 fc 3 + k 2 ) 

1 c 2 c 2 k 2 k 2 (K> + c h k 3 ) (4) _ 2 cy h ktklkl , B) _ 2cik 2 2 kl(3c a h + K>) 



T (3) ± L^^ft-^aV-tv TL ft rv 3 ; (4) z, S^^!^^ T (5) _ 

H 2 K' 2 ' H K' 3 ' X' 4 

2<6) =I W J(7) =I (5) ; (92) 

where AT' := c s fci + Ch(k 2 + fc 3 ). Thus, it turns out that we need to evaluate only Vyjy and Vy^. 

We would now like to define the amplitudes of the above cross bispectra in a similar way as the case of three 
tensors, for which wc have adopted two different normalization conditions, (|74[) and (|88[) . depending on whether 
we are interested in the B-mode polarization or the E-modc polarization and temperature fluctuations. The same 
ambiguity is present for the cases of these cross bispectra, too. Here we simply normalize them in terms of V 2 taking 
into account the fact that these bispectra generate the auto and the cross bispectra of the temperature fluctuation and 
the E-mode polarization, too, which are mainly sourced by the curvature perturbation. Although this normalization 
may not be appropriate for those including the B-mode polarization, we do not touch the issue any further because 
the change of the normalization factor from V 2 to VcVh or V\ can readily be done by multiplying appropriate powers 
of the tensor-to-scalar ratio r. Thus we adopt the following convention: 

R (fW») _ ( 27r ) 4;P C A «hh) , Q „, 

D ij,kl - 1.31.31.3 ^ijM > y^ ) 

where 

7 



4S ] = SH'S^EWiH^i.ka.^'HAi.Aa.fca) + (k a ,M o k 3 ,fc,Z). (94) 

r T c h q=1 

We also define the following cross bispectra: 



(C(ki)f<">(k 2 )f<'»>(k3)) = (2 T r) 3 ( 5(k 1 +k 2 +k 3 ) J B^ ) (k 1 ,k 2 ,k 3 ). (95) 
Here B^ h ^ and are given by 

= -Tjku-A^ = ,3,3,3 T T2r f . 2 E & ^it 3 (ki,k 2 ,k 3 )J^(fc 1 ,fc 2 ,fc 3 ) + (k 2 , S2 o k 3jS3 ), (96) 
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where V^ s3 (ki, k2, k 3 ) is evaluated as 

vg.3 = H - + S3 hf] 2 , v£ 3 = k2 . k3 v2 s3 = fc i- fc |- fc 3 2 v w 3 , 

(3) _ ki-ka (i) _ fef -fc| + fc| (i) (4) _ (i) ( 5 ) _ , 2 V (i) 

^s2,s3 ~~ ^,2 K s2,s3 — 2 |/ s2,s3' K s2,s3 — K s2,s3) "s2,s3 ~ K l ^s2,s3J 

K ( 2 6 ! s3 = ooJL^ i - fc 2 - *s)(*i + fc 2 - fc 3 )(fci - + fc 3 ) [kl - (s 2 k 2 + S3 fc 3 ) 2 ] , V£ s3 = k 2 Mt] s3 . (97) 

C. Two scalars and one tensor 

The cross bispectrum of two scalars and one tensor is given by 

(C(ki)C(k 2 )Mk 3 )) = (2^) 3 5(k! +k 2 +k 3 )B l ( 7 CC ' l) (k 1 ,k 2 ,k 3 ), (98) 



where !?„•£ is of the form 



1 H 6 6 

B ^ h) = 4fc3 fc 3fc3 FlF TC 2 Ch ^cMf{^MM)J {q) {k x MM) + (ki o k 2 ). (99) 



Each contribution is given by 



1 ,.m . .C4i fc 2 



and 



1 cf (fci + fc 2 )(fc 2 + fcifc 2 + fcf) + 2c 2 c h (k 2 + kik 2 + fc|)fc 3 + 2c s c 2 h (k 1 + k 2 )k 2 + c|fc| 



i? 2 if" 2 
1 c 2 fc 2 [2c 2 (fc 2 + gfcifcg + fcf) + 3c a c /l (fc 1 + fc 2 )fc 3 + c 2 fc 2 ] 

1 c 2 c 2 h k 2 k 2 {c s k 1+ K") 

1 c 2 fc 2 [c 2 (fci + fc 2 )(2fci + fe) + 3c s c /l (2fc 1 + fc 2 )fc 3 + 2c 2 k 2 } 



7 (4) 

H A" 3 
J {5) = -Jtj [c 3 s (h + k 2 )(kj + 3fcifc2 + fcf) + 4c 2 c h (fc 2 + 3fcifc 2 + kj)k 3 + 4c s c 2 1 (fc 1 + k 2 )k\ + c 3 h kl] , 
(6) _ 1 cffc 2 fc 2 2 (A" + Qt fc 3 ) 

with AT" := c 8 (ki + k 2 ) + c/jfc 3 . Thus, it turns out that we need to evaluate only . 
As in the case of two tensors and one scalar, we normalize the bispectrum by V 2 as 

,(«/.) _ ( 2 ^) 4 ^C M 

13 h 3 h 3 h 3 V 

where 



B (CCh) _ (2tt) 4 P c 2 (CCfe) 

^i.i ~ ;_3/,3;..3 i.7 ' 1> 1U ^ 



■^f fc) = ^E c X 9) ( k i' k2 ' k 3)^ (9) (fci,fc 2 ,fc 3 ) + (ki ok 2 ), (103) 

We also define the following cross bispectra: 

(C(k 1 )C(k 2 )e (s) (k 3 )) = (2^3<5(k 1 +k2+k 3 ) J B(« /t )(k 1 ,k 2 ,k 3 ). (104) 
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Here B^^^ and A^''' 1 ^ are given by 



= IS? 4CC ' l) = i^^£^:^ c?vi9)(kl ' k2 ' k3)i7(9)(fcl ' fc2 ' fc3) + (kl " k2)> (105) 



where vi 9 ^(ki, k 2 , k 3 ) is evaluated as 



V i X) = S^( fc i - fc 2 - fc 3 )(*i + fc 2 - *3)(fci - fc 2 + fc 3 ), (106) 



and 

V1 2) =V S W, Vi 3 ) = 4v«, VW = SVW, vW=*gvi 1 ) I Vi 6 ) = -|,V«. (107) 



7,2 

Indeed, the above functions are independent of s due to no parity violation. 



1.2" s i "a i»2 8 ' " s "rs i " s 7.2 jL2 



D. Three scalar s 



Here we give the bispectrum defined by 

(C(k 1 )C(k 2 )C(k 3 )) = (2 7 r) 3 5(k 1 + k 2 + k 3 )B^(k u k 2 , k 3 ). (108) 
The result is given in Ref. [H IH[: 



B (CCC) = 

/I iU3 L3 h.3 

Ht/v^ 3 



{k x k 2 k 3 ) 2 C 2 
Cl+ K 



+ c 3 [ E fc' + ^E - ^ E W) + c 4 (e * 3 - 5 E fc ^ 2 - ^ E #3 



9l 

K 2 



2 E k * + E k * k 1 - 3 E rftf - 2fc i fc 2 fc 3 E kik i 

i i=£j i^j i>j 



(109) 



V. EXAMPLES 

In this section, we consider two representative examples of inflation to estimate the amount of non-Gaussianities 
from tensor and scalar perturbations. The first example is general potential-driven inflation studied in Ref. [3{|. This 
class of inflation models includes variants of Higgs inflation enabled by enhancing the effect of Hubble friction. These 
potential driven models have c 2 = 0(1) and c\ ~ 1. Next, to see the impact of generic c 2 more clearly, we study 
k-inflation as another example. 

A. The case of potential-driven inflation models 

We wish to treat a wide class of potential-driven inflation models at one time. For this purpose, we introduce six 
</>-dependent functions to write 

K = -V{4>)+JC{4>)X, G s = h s (ct>)X, Gt = g(<f>) + Iu(<f>)X, G s = h 6 (<f>)X. (110) 

In particular, the above form includes different Higgs inflation models proposed so far (39| . These may also be regarded 
as the Taylor expansion of K(cf>, X) and Gi((f>, X) with respect to X. Would-be leading terms in G 3 and G5 have been 
removed without loss of generality. 
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Slow- roll dynamics of general potential-driven inflation models has been addressed in Ref. [39( . During inflation we 
assume that the following slow-roll conditions are satisfied: 

e=-§,«l, V = ~^«l, S=£- g <l, ^ = m <K ^ (-3,4,5). (Ill) 

It is convenient to define 



u(0:=/C+^, v(<f>):=h 3 + ^, W{ft>):=\ 



Under the slow-roll approximation the gravitational field equations reduce to 

X 



(112) 



6gH 2 ~V, 2e + 5 ~ [u + 3H4>v) . (113) 

Now it is easy to see that Tt — Gt — 2g and 

Ts^^{u + AH^v) ~|(2e + 5)(4-^), Gs ~ -J (u + QH^ru) ~ g(2e + S) (2 - , (114) 

so that J-g and C?s are slow-roll suppressed. It can also be seen that c 2 ~ 1 and c 2 = 0(1). 
The coefficients of the cubic terms are given by 



8iJ 2 24 v ; V VK/ 4 V ; V W 



fc— * ^ 5 ~ ~qM2~-~u/W' E shh ^ w Ch ijE ^ (115) 



and 



Cl * ^(2, + ^) (4- ^ , c 2 ^ -^(2, + *) fl - Ji) + i* fafca - ^ 
1 3 V ; V WJ' 12H y '\ WJ 4H 2 V 65 / 



C5 - Mi 



where /i = ifiXhs as defined in (j35]). It turns out that the other coefficients are of higher order in the slow-roll 
parameter: ~ C3 ~ C4 ~ cq = 0(e 2 ). 

B. The case of k-inflation 

To extract the effect of the nontrivial sound speed, let us consider k-inflation, which is the simplest model with a 
generic value of c 2 . In the case of k-inflation, K = K(<fi,X), G4 = Mpj/2, G3 = = G5, we have 

M 2 e 

T T = Gt = T = M| b rs = M&e, Gs = ~^, 9 = Mp } H, /i = 0, (117) 

c 

with c/j = 1 and r = 16ec s , which simplifies the coefficients in the cubic Lagrangians: 

At 2 f M 2 e 1 

h = b 2 = -^-, bs = --^-, 64 = 65 = 66 = 67 = 0, E shh = —Ch ij E^ j , (118) 

M 2 e 2 M 2 e 2 

Cl =M 2 l£ , c 2 =0, c 3 = ^^, c 4 = c 5 =0, c = -g_ (119) 



^ _ 2H < *" jhi:il ~ 2Hc 2 ^ l ^' j Aa 2 H 2 ^' i ^' j ' 



Note that in deriving the above coefficients we have not invoked the slow-roll expansion. 
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VI. DISCUSSION 



In this paper we have presented the full bispectra, including the cross bispectra of the primordial curvature and 
tensor perturbations, in the generalized G-inflation model which is the most general single-field inflation model with 
the second order equations of motion. 

In the event full observations of these quantities could be made, we could extract many pieces of interesting 
information on the underlying theory. For example, by observing three-point tensor correlation function, we can in 
principle determine the kinetic coupling to the Einstein tensor through /i. Another interesting quantity is the cross 
bispectrum of two tensors and one scalar. If we could observationally identify their coefficients b 2l 63 and b§, we could 
in principle determine Fs, Qs, Ft, and Qt independently with the help of the three-tensor bispectrum which would 
provide a consistency relation of the theory for the tensor-to-scalar ratio ([36| . 

Let us next turn to two-scalar and one-tensor bispectrum whose effective Lagrangian is given by (|53|) . Its most 
interesting component is the first term proportional to c\ = Fs which could be singled out by taking k% small. In the 

standard canonical inflation as well as in k-inflation, the coefficient simply takes ci = Fs = Al^e = -^j — as derived 

in (|119[) . where we have used the consistency relation in the last equality. 

We can also show that this feature remains valid in the case where a sizable local non-Gaussianity is generated as in 
the cases of the curvaton scenario [40l | and the modulated reheating scenarios [ill ] . In such case curvature perturbation 
£ is sourced by another scalar field which we denote by a and its fluctuation by 8a. One can relate £ and 8a as 

C = N a (a)Sa + ^N aa {a)(8a)\ (121) 

using the <5A-formalism (42|. Suppose that a has the Lagrangian C a = K(Y,a) with Y := —(da) 2 /2. Since the 
dynamics of a is practically frozen during inflation and it practically behaves as a massless minimally-coupled field, 
one can expand C a = k(0,<to) + k<t(0,cto)Y in this regime where ao is its expectation value in the domain including 
our horizon today. Then the mean-square fluctuation amplitude of a is given by 

<^) 2 > = a 2 H !n \ = l^ri^' ( 122 ) 
4tt 2 «; ct (0, a ) Nj{a Q ) 

the latter being an outcome of (|121[) , and it determines the relation between 8a and (, too. Then the effective 
Lagrangian representing tcnsor-scalar-scalar coupling is generated from the kinetic term of a in this case and reads 

1 ,„ s,,.„ 1 ( H \ 2 „ Ml x r 



C ssh = -K a (0,a )h^a tfl a tV = - ) T Q h^Q^ u = -^h,,cc,. (123) 

Note that in this case the sound speed is equal to unity. Thus we find that if the sector responsible for the generation 
of curvature perturbations is minimally coupled to gravity with no extra Galileon-like terms, c\ takes the same form 
whether they are generated by the inflaton or another scalar field. Thus this term can provide a test of the generalized 
Galilcon as a source of the structure of the Universe. 

It is a non-trivial issue how to normalize the cross bispectra. In this paper, we have normalized them by the 
power spectrum of the curvature perturbation. This is mainly because these cross bispectra generate the auto- and 
the cross-bispectra of the temperature fluctuation and the E-mode polarization, which are mainly sourced by the 
curvature perturbation. However, such a normalization may be inadequate for the cross bispectra including the B- 
mode polarization. Therefore, we need to directly investigate the impacts on the CMB bispectra [I?], EH . Constraining 
the model parameters by CMB bispectra is a work in progress [H[ . 
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